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Abstract 

We present a new way of quantum kinetic equation derivation. This method appears as a natural 
generahzation of the many-particle quantum hydrodynamic method. Kinetic equations are derived 
,1^ . for different system of particles. First of all we consider quantum plasma and pay special attention 

to the spin evolution. We show that we need a set of two kinetic equations for description of spinning 
particles. One of these equations is the equation for distribution function, however this equation 
contains new function, even in the self-consistent field approximation. This is a spin-distribution 
function introduced in the paper. Therefore we have to derive kinetic equation for spin distribution 
function evolution, which is presented here and used to construct a closed set of kinetic equations. 

^: 

We also present kinetic equation for system of neutral particles with a short-range interaction 
in the first order by the interaction radius approximation. We derive a set of kinetic equations 
for particles having electric dipole moment, this set analogous to the equations set for spinning 
particles, but it has some differences. As a special topic we find kinetic equations for graphene 
carriers in the vicinity of the Dirac points. Derived equations, in general case, contain two-particle 
distribution functions, which take into account contribution of the quantum correlations including 
the exchange interaction, but we restrict ourself by the self-consistent field approximation to obtain 
closed kinetic description, in the system of particles with the short-range interaction. 
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I. INTRODUCTION 



Quantum plasma studying requires developments of new theoretical methods. In classical 
regime it has been very fruitful to use hydrodjTiamical and kinetic methods. So, we can 



expect that their quantum generalization will be as profound. In 



quantum hydrodynamics has been developed by many teams [l| -ll3j |. The semi-relativistic 
interactions, the spin-spin [2|, the spin-current the spin-orbit current-current [l2|] 

interactions, and the Darwin term and the semi-relativistic amendment to the kinetic energy 



ast decade the method of 



12|, have been included in the quantum hydrodynamics scheme along with the Coulomb 



interaction. Wigner kinetic approximation has been used for studying of the quantum plasma 
|l^-j22|. A physical kinetics, free of any special assumptions on the space-time geometry, 
was developed in Ref. 23|]. A quantum distribution function such that calculating statistical 



averages leads to the same local values of the number of particles, the momentum, and the 



energy as those in quantum mechanics was found in Ref. 



24| . This distribution function 



coincides with the Wigner function only for spatia 



24|- 



ly homogeneous systems. It was used 



27j and consideration of some physical 



for derivation of an quantum kinetic equation 

problems 25|- 29j. Kinetic equation for charged spinning particles defining the distribution 



function of ten variables the coordinate r, the momentum p, the spin s, and the time t 
was considered in Ref. jl4|. Integration of this distribution function /(r, p,s,t) on the spin 
and momentum gives concentration of particles n(r,t). Integrating the product of the spin 
vector and the distribution function on the same variables gives spin density. Quantum 
kinetic equation forf (r, p, s, t) with the semi-relativistic effects caused by interaction was 



suggested in Ref. 



181]. Obtaining of the kinetic equation by averaging of the one-particle 



hydrodynamic equations on the Maxwell or the Fermi distribution function was considered 
in Ref. [s^. Set of equations consisting of the Vlasov equation, including ponderomotive 
force acting of the magnetic moments, for particle motion and Bargman-Michel-Telegdi 
(BMT) equation, describing magnetic moment evolution, was used in Ref. 3l| to consider 
dispersion of elementary excitations for the tensor gyromagnetic ratio. Review of the recent 
achievements in the quantum kinetics can be found in Ref.s 32|, 33 1. 

Microscopic density of particles in classical physics can be presented as the sum of Dirac 



2 



delta functions j34l|-|36| 



N 



r,t) = 5^5(r-r„(t)), 



n=l 



where is the total number of particles in the system. This definition of the particles con- 
centration was used for construction of the quantum particles concentration and derivation 
of the quantum hydro dynamic equations l|,]8i, jsl- Corresponding microscopic distribution 
function can be written as 



/ = 5^5(r-r„(t))5(p-p„,(t)). 



Averaging of this function leads to macroscopic distribution function allowing to derive a 



chain of relativistic kinetic equations 



General form of quantum distribution function 



might be presented in the following form 2A\ 

/(r,p) = ^nn{v)nn{r), 

n 

and it was expressed via the one-particle statistical operator. However, we in this paper 



consider the quantum mechanical averaging 



<L>-- 



^p*L^pdR 



of / defined as 



f = ^Sir -rn)S{p-pJ, 



(1) 



where p„ is the momentum operator for n-th particle, and p is the numerical vector function 
which arithmetizes the momentum space, as the coordinate r arithmetizes coordinate space. 

We show that evolution of the quantum mechanical averaging of the operator ([1]) leads to 
arising of the second function, which is the quantum mechanical averaging of the following 
operator 

^" = 5^a:5(r-r„)5(p-pJ. (2) 

n 

Thus, we should use a set of two kinetic equations for description of spinning particles 
evolution. 

We use this method for derivation of kinetic equations for different physical systems. 
These are the spinning quantum plasma briefly discussed above, charged or neutral particles 
having electric dipole moment, neutral particles with the short-range interaction, graphene 



carriers and graphene excitons. Corresponding quantum hydrodynamic equation have been 
obtained previously. At both the graphene carriers and the graphene excitons description we 
have deal with the quasi-spin. Consequently, we have to obtain a couple kinetic equations 
as for spinning quantum plasma. At description of the electrically polarized particles lead to 
appearing of a new function, analogous to ([2]), but we should put the electric dipole moment 
operator instead of the spin operator. 

This paper is organized as follows. In Sec. II we start our microscopic derivation of 
quantum kinetic equation. This derivation is presented for the simples case of charged par- 
ticles in an external electric field described by scalar potential with the Coulomb interaction. 
The self-consistent field approximation considered for obtained kinetic equation. In Sec. Ill 
we introduce the distribution function in a magnetic field and derive kinetic equation for 
charged spinless particles in an external electromagnetic field. In Sec. IV we consider kinet- 
ics of spinning particles, we derive a couple of kinetic equations for the distribution function 
and spin distribution function. In Sec. V we obtain kinetic equations for particles having 
electric dipole moment. In Sec. VI we study kinetic of neutral particles with the short range 
interaction. In Sec. VII we find a set of kinetic equation for carriers in the graphene. In 
Sec. VIII a set of kinetic equations for the graphene excitons is obtained. In Sec. IX we 
present the brief summary of our results. In Appendix we critically examined some general 
method derivation of kinetic equation used as for quantum and for classic systems. 

II. CONSTRUCTION OF MACROSCOPIC EQUATIONS 
A. Kinetic equation for spinless particles: General form 

The equation of quantum kinetics is derived from the non-stationary Schrodinger equation 
for system of N particles: 



The following designations are used in the equation ([3]): = —ihd" — tn-^^.extl ^n,ext, 
gj,^ are the potentials of the external electromagnetic field, (9° = is the derivative on 





(3) 
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the space variables of n-th particle, and Gnk = '^/fnk is the Green functions of the Coulomb 
interaction, tp{R,t) is the psi function of N particle system, R = (ri, rjv), e„, m„ are the 
charge and the mass of particle, h is the Planck constant and c is the speed of light. 

The first step in derivation of the kinetic equation is a definition of distribution function 
of particles. We determine the distribution function of particles as the quantum-mechanical 
average of the operator of distribution function: 

/ = -rn)S{p-Pn)- 

n 

This function is the microscopic distribution function in classic physics, 5(r)-is the Dirac's 
(5-function. 

In that way, the distribution function has form: 

/(r,p,t) = ^ I ^*(i?,t)5^('5(r-r„)(5(p-pJ + 5(p-pJ5(r-r„,)')v;(i?,t)rfi?, (4) 

where dR = Y[n=i ^^n- This definition symmetric relatively to operators S{r — r,„) and 
6{p — Pn), but it still is not fully symmetric. To get the fully symmetric definition of the 
distribution function we need to add complex conjugated quantity. Thus we have following 
definition 

/(r,p,t) = ^J (^r{R,t)Y,(^^i^-^n)S{p-pJ+S{p-pM^-r^^^ (5) 

where c.c. stands for the complex conjugation. 

For the first step we will consider set of the charged particles in external quasi-static 
electric field including the Coulomb interaction between particles, this means that we do 
not include the vector potential. Next, we will consider electromagnetic external field, and 
we will also consider spinning particles with spin-spin interaction. 

Integrating of the distribution function over momentum we have concentration 



n{r,t) = J /(r,p,t)rfp (6) 
where the concentration has familiar from the quantum hydrodynamics form [sl, [l| 

n{r,t) = j dRY,K^-rn)r{R,t)tP{R,t), (7) 
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Integrating of the product of momentum on the distribution function over momentum we 
find particles current 

j(r,t) = J p/(r,p,t)rfp, (8) 

details of calculations you can find in appendix. 
In the same way we have the kinetic energy 

f 

£fcm(r,t)= / — /(r,p,t)rfp. (9) 
In the absence of the inter-particle interaction we find 

dtf + -pdrf + e^^ir, t) sin(^rVp)/ = 0, (10) 



where 

sm(^,Vp) = f |5£^(^,Vp)-«. 

In equation ( fTOl) we have used designation (p{r,t)'^r, it means that the gradient operatorV 
on spatial variables acts on the left-hand side, instead of usual acting of operators on function 
standing on the right-hand side. We do not write the Plank constant h in the argument of 
sin to make this notation more handy. 

The last term appears here due to the commutation of the electric scalar potential with 
the 5(p — p„). At / > 2 terms are proportional to the plank constant ~ h''~^. They disappear 
in the classical limit h ^ 0. 

In the quasi-classical limit we have to include one term of the sum in the last term only, 
so we have 

dtf + -pdj + eEVp/ = 0, (11) 

m 

where we have written E = —dr^p. 

We also admit that distribution function (j4]) satisfy to the following equation 

dtf + -pdrf + e'-ifiv, t) J2 ^(VrVp)7 = 0. 

1=1 

At derivation of kinetic equation the distribution current arises in the second terms of 
equations ffTOj) and f fTOjl 

J"(r, p, t) = 1 y f {p^riR, t)) (^(r - rn)5(p - pj + S{p - pJ5(r - r„)) ^R, t) 



+ r {R, t) (^(r - r„)5(p - p J + 5(p - p J5(r - v,S^py{R, t) + c.c. j dR, (12) 

which we approximately present in the following 

J°(r,p,t)=p"/(r,p,t), (13) 

which corresponds to the traditional structure of kinetic equation. The last formula has 
been already used in formulas ffTOj) and f|T9l) . 

Integrating equation (fTO!) and other form of kinetic equation obtained in this paper we can 
find that they give us the equations of the many-particle quantum hydrodynamics obtained 
in Ref.s Q, fl, 

Now we are going to include the Coulomb inter-particle interaction 

dtf + — P^r/ + e^v? sin(^rVp)/ 

X (^r {R. t) (r - r„)5(p - p J + <5(p - p J(5(r - r„)^ i^{R, t) + c.c. j rfi? = 0. (14) 

In equation f[T^ we find one more term in comparison with equation fITU]) . This new term 
is the last term in equation (fT^ . it is caused by the Coulomb interaction. 
Equation (fT4|) can be rewritten using two-particle function 

•9*/ + — P<9r/ + e^v^ sin(^rVp)/ 
m ft 

+ e^^ y G(r - r') sin(^rVp)/2(r, p, r', p', t)dv'd^' = 0, (15) 

/2(r,p,r',p',t) = l jiriR.t) 5(r'-r,)5(p'-p,)x 

X (^(5(r - r„)(5(p - pj + (5(p - pj(5(r - r„)^ ^/;(i?, t) + c.c. j di?, (16) 

is the two-particle distribution function. 

In the quasi-classical limit we have to include one term of the sum in the last term only, 
so we have 

dtf + -pdrf + eEVp/ - / VrG{r - r') ■ V^r, p, r', p', t)dr'dp' = 0. (17) 
m J 



where 



Equation (fTTI) looks like classical kinetic equation. However we should admit that two- 
particle distribution function /2 is defined via wave function ip. Therefore, we see, analo- 
gously to the quantum hydrodynamics jsl, this equation is a quantum kinetic equation 
containing information about quantum effects, in particular, the exchange interaction. Some 



methods o: 
in Ref. 3^ 



obtaining of the closed quantum kinetic description was discussed and developed 



B. The self-consistent field approximation 

Considering many-particle quantum kinetics we find that equation of distribution function 
/ evolution contains two-particle distribution function /2 that corresponds to the classical 
kinetics. For particles with a long-range interaction, in our case it is the Coulomb interaction, 
we can introduce, well-known in classical physics and quantum hydrodynamics, the self- 
consistent field approximation. In this approximation we separate two-particle distribution 
function /2 in the product of two one-particle distribution functions / 

/2(r, p, r', p', t) = /(r, p, t)f{r', p' , t). (18) 

This approximation gives us closed mathematical apparatus in the form of one nonlinear 
integral equation, which in the quasi-classic limit to be 

dtf + -pdj + eEVp/ - / " r')/(r , P , t)dr'dp' ■ Vp/(r, p, t) = 0. (19) 

m J 

In general case, from equation (120|1 . in the self-consistent field approximation, we find 

dtf + —pdrf + e^ip sin(^rVp)/ 
m n 

+ e^l J G(r - r') sin(^Vp)/(r, p, t)/(r', p', t)dr'dp' = 0, (20) 

Introducing electric field for inter-particle interaction we come to the set of kinetic equa- 
tion and field equations (Maxwell equations) instead of integro-differential equation (IT^ . 
This set appears as 

dtf + -pdj + eEVp/ = 0, (21) 

m 

VxE = 0, VE = 47re//(r,p,t)dp. (22) 

In the result we see that derived equations corresponds to the Vlasov equation. It also 
coincides with the results obtained in Ref.s 15|- 2l|, where was used the Wigner distribution 
function. 
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III. CHARGED SPINLESS PARTICLES IN MAGNETIC FIELD 



Studying charged particles in an magnetic field we have long derivative D„ = —iKWn — 
eKnjc instead of the short one —ihSJ^i- So, we have to include it at definition of the 
distribution function /(r, p, t) and in this case we have 

/(r, P,t) = \j (^i^* t) (^(^ - ^-)'^(P - + - ^-)'^(^ - ^")) ^) + ^-^^ j 

(23) 

In the absence of the inter-particle interaction we find 

dtf + -Pdj - —d^A^d^^ (py\ - -(9*A)Vp/ 
m mc \ J c 

- -^lA-sin(^rVp)(p"/) + e^V'sin(^rVp)/ = 0. (24) 

Including in consideration of the vector potential of the external electromagnetic field leads 
to appearing of three additional terms in the kinetic equation (!24|) . which are the third- fifth 
terms. The fourth term gives contribution in the force acting on the charge of the external 
electric field. The third and fifth terms is the magnetic part of the quantum Lorentz force. 
In the quasi-classical limit it has well-known form presented below. 

In the quasi-classical limit we have to include one term of the sum in the last term only, 
so we have 

dtf + -drf + e f E + — p X Vp/ = 0, (25) 
m \ mc J 

where we have written E = — SrV? — dtA./c. 

We should write the Lorentz force in equation ( 125|) via velocity p/m = v, but distribution 
function / = /(r, p, t) depends on momentum, and we write all coefficients via momentum 
P 



IV. KINETIC EVOLUTION OF SPINNING PARTICLES 

In section III we present the Schrodinger equation contained only the Coulomb interac- 
tion. Now we going to consider kinetics of spinning charged particles. In this section we are 
particularly interested in the spin-spin interactions. The whole Hamiltonian, used in this 
case reads 



9 



\ i^^^nGun - l^luGtalai), (26) 



2 



where 



In these formulas we have used following designations. The Green's functions of the 
Coulomb, the spin-spin and the spin-current interactions has the following form Gpn = l/vpn, 
Gpn = 47r5°'^5(rp„) + 9p (9^(l/rp„), where 7p is the gyromagnetic ratio. For electrons 7p reads 
Ip = Cph / {2mpc) , Cp = — |e|. The quantities ipp^^ = (p{rp,t), A'^f^^^^-^ = A°'{rp,t) are the scalar 
and the vector potentials of the external electromagnetic field: 



B^e.t)M = e-^^V^pAl^,^{rp,t), 
Ete.t)M = -V^^e.tM - ldtA:^,irp,t). 



Gp is the Pauli matrix, a commutation relations for them is 



\a^^]=2i5p^e''P^a;. 



First we present kinetic equation for spinning particles in the absence of the inter-particle 
interaction 

dtf + -Vdrf - —d^A^dl fp"/) - -(9*A)Vp/ 
m mc \ / c 

- — ^A"sin(^rVp)(p°/) + e^y.sin(^rVp)/ 

-/i^5"sin(^rVp)5"(r,p,t) = (27) 

In the quasi-classical limit we have to include one term of the sum in the last term only, so 
we have 



dtf + -dj + e f E + — p X Vp/ + d^B^{r, t)dp^S^{r, p, t) = 0, 
m \ mc J 



(28) 



where we have new quantity in these equations, we can call it spin-distribution function, it's 
explicit form is 

5"(r,p,t) = (^r{R,t)Y,(^^i^-^n)S{p-B^)+S{p-BMr-rn)^ 

(29) 
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where h.c. stands for the Hermitian conjugation. 5*' 
density, which arises in the quantum hydrodynamics |2| 



t) is kinetic analog of the spin 
^ and has form 



(30) 



we have used same letter for designation of the spin density and the spin-distribution func- 
tion, but they differ by set of arguments. 
It can be shown that 

^"(r,t) = J ^-(r,p,t)rfp. (31) 
Magnetization M"(r,t) usually used in the quantum hydrodjTiamics j^, [J], jo], and 60 |. 



Magnetization M°'(r,t) has simple connection with the spin density M°(r, t) = 'jS°'{r,t), 
where 7 is the gyromagnetic ratio for considering species of particles. 

To get complete description of systems of spinning particles we have to derive an equation 
for the spin-distribution function S"'{r,p,t). 



A. Including of interaction 

Introducing of interaction makes the kinetic equation more larger, but we present it's 
explicit form 

dtf + -pdj - —d'^A^d^ (p'^f) - -(9iA)Vp/ 
m mc ^ \ J c 

- — ^A"sin(^rVp)(p"/) + e^<^sin(VrVp)/ 
mc n n 



-/i^5"sin(^rVp)5"(r,p,t) 

^\'n j jenCkGnk) sin(^rnVp) 



X ( riR, t) ( ^(r - r„)5(p - D„) + 5(p - D„)5(r - r„)J ^P{R, t) + h.c. ) dR 

n,k^n 

X I r{R. t) ( 5(r - rn)5(p - D„) + 5(p - D„)5(r - vS] «V'(^, t) + h.c. ]dR = Q. (32) 
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Let's describe meaning of terms in this equation. The first seven terms coincide with the 
same terms in the kinetic equation for charged spinless particles i^^. They are time evo- 
lution of distribution function presented by the first term, convective part of distribution 
function evolution, presented by the second and third terms, correspondingly. Terms four, 
and six give quantum-kinetic generalization of the Lorentz force describing action of the 
external magnetic field on particle charges. The fifth and seventh terms show the force 
acting on charges exerted by the external electric field. The eighth term describes action of 
the external magnetic field on spin of particles entering equation via the spin-distribution 
function S'"(r, p, t). The last two terms present inter-particle interaction, the Coulomb and 
the spin-spin interactions correspondingly, and it's influence on evolution of the distribution 
function. 

Introducing two-particle spin-distribution function as 



^r(r,p,r',p',t) = i jir{R,t) 5(r'-r,,)5(p'-D,)x 



X (^5(r - r05(p - D„) + 5(p - D„)5(r - r„) J t) + h.c. ) dR (33) 

we can rewrite previous equation in the following form 

dtf + -Pdj - —d^A^d^^ fp"/) - -(5*A)Vp/ 
m mc ^ \ J c 

- — ^^"sin(VrVp)(p"/) + e^y.sin(VrVp)/ 
mc n n 

-/i^fi'^sin(^rVp)5"(r,p,t) 
+e^'- j G(r-r')sin(^rVp)/2(r,p,r',p',t)rfr'rfp' 

+ 7'^ j G^''(t - v') sin(^rVp)Sr(r, p, r', p', t)dv'dp' = (34) 

Using of the two-particle functions we have the two last terms in equation (IMj) written 
shortly. It also gives a form of kinetic equation analogous to the BBGKY equations, that 
allows us use well-known and useful ideas to get closed set of equations, such as introduction 
of the self-consistent field approximation, which we discuss for spinning below. 
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B. spin-distribution function evolution 



In kinetic equation for spinning particles appears the spin distribution function. There- 
fore, for construction of the closed set of equation describing spinning particles we have 
to find equation evolution of the spin distribution function. For this goal we differentiate 
spin distribution function with respect to time, after some calculations we find the kinetic 
equation for spin distribution function evolution 

5t5"(r,p,t) + — p9r5" - —d"'A^d^(p"'sA - -(9tA)Vp5" 
m mc ^ \ J c 

- — 4A^sin(^rVp)(/S") + e4</^sin(Wp)^" 
mc n n 

-/x^E°sin(Wp)/(r,p,t) 
+e^^ j G{r - r') sin(^Vp)M2(r, p, r', p', t)dr'dp' 
I G"^(r-r') sin(^rVp)iV2^(r,p,r',p',t)dr'dp' 

_ ^ ^ I G/^^(r,r')S2^'(r,p,r',p',t)dr'dp'j. (35) 

There is no high order derivatives on space variables in the last term. In the second term of 
equation fl35l) we have used following presentation for the kinetic spin current 

J°^(r,p,t)=/5"(r,p,t), (36) 

which has following explicit form 

J-^(r, P,t) = ^J (^{D^^riR, t)) (^(r - rn)5(p - + 5(p - D„)5(r - r„)^ ^^^(i?, t) 

+ r {R, t) (^(r - r„)5(p - D J + 5(p - D J5(r - r„)^ D^a^iR, t) + /i.e. j dR. (37) 
We have used a couple of new two-particle functions 

M2"(r,p,r',p',t) = iy"u*(i?,t) 5(r'-r,.)5(p'-D,)x 
X f 5(r - r„)5(p - D„) + 5(p - D„)5(r - r„)^ a'^^{R, t) + h.c. ] dR, (38) 
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and 

N^{T,p,r',p',t) = ^ j(r{R,t) J2 S{r'-n)S{p'-f>k)x 

X (^{r - r„)5(p - D„) + 5(p - D„)5(r - r„)^ t) + /i.e. j dR. (39) 

These two two-particle distribution function describe correlations between spin dynamic of 
one particle with the space motion (concentration evolution) other particle. Consequently, 
in the self-consistent field approximation this function represents as 

M^{v, p, r', p', t) = 5°(r, p, t)/(r', p', t), (40) 

and 

iV2"(r,p,r',p',t) = /(r,p,t)5"(r',p',t). (41) 

Using equation fl35|l we can find hydrodynamical equation for spin evolution S'"(r, t) = 
^"(r, p, t). Terms containing Vp/(r,p,t) give no contribution in the hydrodynamical spin 
evolution equation. Consequently, the First, second, and last terms of equation fl35l) give 
contribution in the hydrodynamical spin evolution equation, and we have 

9i5"(r,t) + - / ar(p5"(r,p,t))dp 
m J 

- '^e-P^ (^^S\r, t)+^ j G^\r, r')Sf{r, r', t)rfr' j = 0, (42) 

which corresponds to equations presented in Ref.s ^, 9|. We admit that we have not 
introduced velocity field in equation P2|) . The second term in equation fH2|) is the divergence 
of the spin current and the last term describes interactions. Hydrodynamical definition of 
spin-current we briefly discus in the Appendix, where we show that presented here definition 
corresponds to the hydrodynamical one. 

The second term in equation fH21) is the divergence of the spin-current, and the last term 
in a big brackets describes interaction. Hydrodynamical spin-current definition we briefly 
discuss in Appendix. 



C. selfconsistent field approximation for spinning particles 

In previous sections we have considered the self-consistent field approximation for charged 
spinless particles. After derivation of the general set of kinetic equations for spinning par- 
ticles we should consider self-consistent field approximation of these equations, to be short, 
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we present equations in the quasi-classic limit. We have presented approximate form for 
and A'^^, but we have one more two-particle function , which is represented as 

Sf{r,p,r',p',t) = S^{T,p,t)S^{T',p',t) (43) 

In the result we have next set of equations 

dtf + -dj + e f E + — p X Vp/ + d^B^{r, t)9p„^^(r, p, t) = 0, (44) 
m \ mc J 

and 

9*5" + — 9r^" + e( E + — pxB )Vp'5" + 9^5"(r,t)9p^/(r,p,t)-^5"^^5^^^ = 0, (45) 
m \ mc J 

Due to the fact that we have considered spin-spin interaction and we have not included spin- 
current and current-current interaction we have magnetic field satisfying to the following 
equation V x B = 47rV x J S{r,p,t)dp. Using additivity of the electromagnetic fields we 
can include magnetic field caused by the electric currents, therefore we have 

V X E = 0, VE = 47re / /(r, p, t)dp (46) 

V X B = 47rV X / S(r, p, t)dp + / p/(r, p, t)dp, VB = 0. (47) 

Equations dH]) and ( 145|) are the set of kinetic equations for the quantum plasma of 
spinning particles. Kinetic equations for quantum spin plasma have been derived in other 
papers 15|- 22| However, obtained in this paper equations have several differences with 
attained before, and main difference is the fact that we get set of two kinetic equation for 
distribution function / and the spin distribution function S*", instead of one equation for 
the distribution function / only. 



V. KINETICS OF PARTICLES HAVING ELECTRIC DIPOLE MOMENT 

The electric dipole moment is the property of neutral and charged particles which is, 
analogously to the spin, leads to anisotropic inter-particles interaction. Analogously to the 
spin kinetics we have to have two kinetic equations, one for usual distribution function (|5]), 
and another one we can called the dipole distribution function, which is the average of the 
following operator 

P" = ^C5(r-?„)5(p-pJ, (48) 

n 
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where dn is the electric dipole moment of particle, usually molecules happen to be polar 
particle, i.e. having electric dipole moment. Quantum hydrodynamics of particles having 
electric dipole moment was considered in Ref. jsl. Boltzmann equation has been used 



in Ref.s 



40|. As an example of infinite bubble diagram expansion method using for 



evolution of particles having electric dipole moment see Ref. [4l|. To reveal main features 
of the electric dipole moment kinetics we present kinetic equations for neutral particles 
having electric dipole moment. We start our derivation from the many-particle Schrodinger 
equation with the Hamiltonian 



H 



E (^P? - d-El.t + ^t.ap(r., t) j - \ g [d-d'.C:^^ , (49) 



where the Hamiltonian of the electric dipole interaction is 

,1 ..^p 



Hm = -d^d"- ■ d'^dl (50) 



which corresponds to the Maxwell's equation, as it was shown in Ref. 42 1. 
Using well-known identity 

_ q^qpI ^ - 3ryy ^ ^s-P6{t), (51) 
r r"^ 3 

we can see that the Hamiltonian f l50|) differs from usually used one 

H^. '"-';:^'/^' ,y^,. (52) 

Necessity to consider the Hamiltonian of electric dipole interaction in the form fl50l) caused 
by the fact that it must accord to the Maxwell's equation. This connection exists since 
Maxwell's equations describe electric field and it's connection with the sources. In our case 
source is the density of electric polarization. Action of the electric field on the density of 
polarization come in to equations of motion via the force field. As the result it describes 
interaction of polarization (electric dipole moments). To be short we present them in the 
self-consistent field approximation 



and 



dtf + -drf + e f E + — p X Vp/ + 9,5^(r, t)9p,P^(r, p, t) = 0, 
m \ mc J 

dtP'' + -drP"" + e ( E + — p X B ) VpP" + d^E^ir, t)9p-,P"^(r, p, t) = 0, 
m \ mc J 



(53) 



(54) 
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where 

R-^ir, P,t) = \j E Cc?^(5(r - ?n)5(p - P J + 5(P - P J5(r + c.c. \ dR 

is the one particle distribution function depending of product of the two operators of the 
electric dipole moment of one particle, and the electric field caused inter-particle interaction 
satisfy to the Maxwell equations 

V X E = 0, VE = -47rV/P(r,p,t)rfp. (55) 

At development of the quantum hydrodynamical description of systems of particles having 
electric dipole moment js] has been suggested a method of closing of the QHD set. Analogous 
approximation can be used here. 



VI. KINETICS OF NEUTRAL PARTICLES WITH SHORT RANGE INTERAC- 
TION 

Presented derivation can be used for systems of particles with the short-range interaction. 
The quantum hydro dynamics of neutral particles with the short-range interaction was con- 



sidered in Ref. [43|, where ultracold quantum gases were considered particularly. Following 



ideas developed in this paper and in Ref. 43| we can derive kinetic equation for described 



systems. To be short we present corresponding kinetic equation in the quasi-classical limit 

dtf + -dj - VV^^tVpf - [ Vrf/(r - r') ■ Vp/2(r, p, r', p', t)dr'dp' = 0, (56) 
m J 

using the fact that we consider the short-range interaction we can rewrite general equation 
( 156|1 in the following form 

dtf + -dj - VKxtVp/ + Ird^dp [ Ur, p, r, p', t)dp' = 0, (57) 
m 2 J 



where 



r.ff,rirf-m, (58) 



and we have made additional approximation, we have considered the short-range interaction 



in the first order by the interaction radius 43|. In general case the short-range interac 



tion leads to the expansion of the term describing interaction in the series, and we have 
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included only first term of this expansion. This scheme was developed in quantum hydro- 



dynamics 



43|, contribution of the next nonzero term in this expansion, in hydrodynamical 



approximation, can be found in Ref.s 



43j, |4J], |45|, |46|. 



Integrating the last term in equation flFTI) gives us trace of hydrodynamical two-particle 
concentration at great length considered in Ref.s 43|, 44 1 



n2(r, r, t) = Trn2(r, r, = J /2(r, p, r, p, t)dpdp'. 



Trace of the two-particle concentration was calculated in Ref. 43|], where was shown, in 
the first order by the interaction radius, that 

n2(r,r,t) = n^(r,t), 

for system of particles being in the Bose-Einstein condensate state. 



Methods developed in Ref. 



43| can be used for further studying of the structure of the 



quantum kinetic equation for quantum gases. 



VII. GRAPHENE KINETICS 



This section is dedicated to development of the quantum kinetic theory for graphene 
carriers. It is important to mention several interesting papers have been recently published, 
where possible applications of graphene for future devises is discussed |47t|, [48| and history 



of graphene studies is presented 

Methods of physical kinetics are usually used for studying of collective properties of 
many-particle systems, such as dispersion dependencies of collective excitations and their 
linear and nonlinear evolution. Properties of collective excitations in graphene, including 



50|, [5l|, |52|, |53| by means of the 



spectrum of linear plasmons, were considered in Ref.s 
random phase approximation (RPA). 

Working with the graphene we have deal with the system of charged spinning particles 
in the presence, in general case, of an external magnetic field. Therefore we have to use the 
distribution function defined by formula (!23|) . Before derivation of the kinetic equation for 
graphene we will describe basic equation which drives graphene electrons in the vicinity of 
the Dirac points. 
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We use the many-particle spinor massless Dirac equation [5J], [55 1 

ihdtil) = (^vpa^Di + Cicpi^ext^ + ^ ^ eiCj-Gij^ tp. (59) 

The following designations are used in the Hamiltonian flM?]) : -D" = —ihd'j' — eiA^^^^/c, (pi^ext, 
^fext ~ is potentials of the external electromagnetic field, Ej e^.^ = —Vipi^ext — dt-^i,ext is 
the electric field, 'Qi^ext = curlAi ext is the magnetic field, quantities Cj, mj-are the charge 
and mass of particles, h-is the Planck constant, and Gij = l/vij is the Green functions of 
the Coulomb interaction. In equation ( 159|) the spinor wave function ip = ilj{R,t) depend on 
2N coordinates R = [ri,...,rAr] and time, where rj = is the 2D coordinates of each 

particle. Potentials (fii^ext = fext{''^i,t), ^text — '^exti^iy'^) ^i^o depend on 2D variables. This 
fact has deep consequences. Potential part of electric field connected with the scalar potential 
via space derivative: Ej = — VjV^j. Consequently in equation f l59|) there is no contribution 
of external electric field directed perpendicular to the graphene plane (Contribution of 
Ez might appear via dtAz). Physically, there is no limitation on attendance of z projection 
of electric field and it's action on graphene electrons. The magnetic field vector to be 

B = curl A = ex{dyAz — dzAy) 

+ey{dzA^ - d^Az) + ez{d^Ay - dyA^), 

two component of the vector potential of the magnetic field A^., Ayaie presented in Hamil- 
tonian fl59|) . and they does not depend on coordinate z. Therefore equation fl59|) contain z 
component of the magnetic field only. In this paper we interested in action of the exter- 
nal magnetic field directed at angle of graphene plane. Therefore, we get graphene kinetic 
equations including whole vector of magnetic field i? ■ — ?■ B = [i?^-, By^ Bz]. 

Quantum hydrodynamics corresponding to described microscopic theory was developed 



in Ref . 



58 



for derivation of the 



561], [57|]. Kinetic equation for graphene was used in Ref. 
hydrodynamic equations. 

It should be admitted that the spin matrix vector a" describe not spin, but so called quasi- 
spin. This notion appears in graphene due to the fact that that Wigner-Seitz cell lattice 
contains two carbon atoms. Thus, graphene lattice consists of two sub-lattices, which are 
usually called lattice A and lattice B. Therefore, in his motion the graphene pz electron 
moves from atom A (B) to atom B (A), and transitions of the graphene electron between 
sub-lattices A and B is modelled by means quasi-spin. 
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Kinematic term in equation ( 159|) is the helicity of particles a^p", so we have here that the 
Hamiltonian is proportional to the helicity of particles, and we consider dynamics governed 
by the spin 1/2 two dimensional helicity as the kinetic energy of particle. In Ref. 



we 



have considered the QHD formulation of the particle dynamics governed by the spin-1 two 
dimensional helicity. It has been suggested that such equations might describe evolution of 
excitons in the graphene, existing in two layer graphene. In the next section we will consider 
kinetic formulation of this model. Here we continue development of the kinetic theory for 
electron in monolayer graphene. 

Differentiating distribution function ( 123|) with respect to time and using equation ( 159|) 
for time derivatives of the wave function we find the graphene kinetic equation. In general 
case it has form 



-^F-J^A^ sin(^rVp)5" + e^(^ sin(^rVp)/ 

+ e^'- j G(r,r')sin(^rVp)/2(r,p,r',p',t)(ir'(ip' = 0. (60) 

This equation contains the spin distribution function ^"(r, p, t), so we will derive equation for 
the spin distribution function S'°(r, p, t). As well, we rewrite equation (162!) in the classical- 
like limit 

dtf + vfVS + eEVp/ - vp-e'^^^B^dl^S'^ 



- e'f ^"^(r, r')<9^/2(r, p, r', p', t)dr'dp' = 0. (61) 



As it discussed in Ref. 



50|, graphene properties have no classic limit, and the notion 



"classical- like limit" is used here to show that we consider the first term of operator 
sin(^rVp) only. Obtained equation flMjl is quantum equation, but it accounts leading 
on h terms only. 

We have to admit that vectors r and p are two-dimensional vectors, however vectors 
E(r, t) and B(r, t), and even pseudo-spin vector S(r, t) can be considered as three dimen- 
sional vectors. Belonging of the graphene electrons to the lattice A and B is modelled by 
increasing of amplitude of tpA or tpB, where tp* = {ipA^i^B) spinor function of electron, and 
appearing of the Sz during derivative at using of the commutation relation for spin matrices 
is quite reasonable. Since belonging of the electron to lattice A or B describes by the struc- 
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ture of the Pauli matrices, and we can easily use some combination of the Pauh matrices as 
they appear at derivation. 

Kinetic equation for the spin distribution function appears as 

Ti c ^ 

-VF-d^A^d^f + vf-Ia'' sin(^rVp)/ + e^cp sin(^rVp)5" 

+ e^^ J G(r, r') sin(^rVp)M2"(r, p, r', p', t)dr'dp' = 0, (62) 
and in the classical like limit it has form 

dtS" + VF&^f + '^e'^^^p^S^ + eE^d^S'^ - VF-e"^^BW^f 
n ^ c 

-e^ J VrG(r, r') VpM2"(r, p, r', p', t)dr'dp' = 0. (63) 

In the absence of inter-particle interaction we can rewrite the couple of graphene kinetic 
equations 

dtf + vfVS + eEVp/ - VF-e^'^^B^d^S'' = 0, (64) 



e 
c 

and 



dtS"" + VFd^f + '^e'^P^p^S^ + eEVp^" - VF^e'^^^BW^f = 0. (65) 

In the self-consistent field approximation the graphene kinetic equation (!6T]) . ( 163|) trans- 
forms in 

dtf + VfVS + eEVp/ - VF-e"^''B^d^S^ 



e 
c 

2 



e 



J 9"G(r, r')/(r', p', t)dr'dp' ■ 9^/(r, p, t) = 0, (66) 



and 



e 



2vi 

2 



j VrG(r, r')/(r', p', t)dT'dp' ■ Vp^"(r, p, t) = 0. (67) 

Due to two-dimensionality of the system we can not introduce the field of interaction, or if 
we actually would be forced to do it we get S function in the right-hand side of the field 
equation VE= AiTp= ATre6{z) J /(r, p, t)(i^p. 

In the result we have the closed set of quantum kinetic equations for graphene, where 
evolution of the two distribution functions obeys to the two kinetic equations, see for example 
dMland dSZ])). 
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VIII. GRAPHENE EXCITON KINETICS 



In the previous section we have considered kinetic equation for carriers in graphene, 
which can be considered as a spin-1/2 hehcity governed system of charged particles. In this 
section we going to study kinetic properties of the spin-1 helicity governed particles motion. 
This model has been suggested in Ref. 59| for obtaining of the quantum hydro dynamical 
descriptions of the graphene excitons and dispersion dependence of collective excitations 
of the Bose-Einstein condensate of the graphene excitons. Here we present corresponding 
kinetic equations. 

We present a basic equation here for evolution description of excitons in graphene 



(68) 



Using this equation we derive and present below the equations for description of excitons 
collective motion. Equation fl^S]) differs from flS9]) by the form of spin matrix and the form 
of interaction. Equation f l68|) contains following quantities: wave function ip = ip{R,t), R is 
the whole particles coordinates R = [ri, r^, ttv], = [xj, yj, sf are the spin-1 matrixes 
for i-th particle, = —iKW is the momentum operator, Vi^ext is the potential of external 
field, Uij is the short-range interaction potential describing the interaction between excitons 
in graphene. We notice that at the same time some particle might interact with several 
particles, by means short-range interaction potential Uij. For this statement illustration 
we refer to the liquid where molecules is neutral and interacts with the several neighbor 
molecules. We consider quasi spin-1 particles and spin operators are 3x3 matrixes 



'^O 1 0^ 



1 

V2 



1 1 
1 



1 

V2 



( 



\ 



-z 

1 -I 
% 



1 ^ 


-1 y 



and the commutation relation for spin-1 matrixes is 



(69) 



Differentiating the distribution function , defined by formula ( !23l) . but without vector 
potential, with respect to time we find following kinetic equation 



dtj + vpVS - VKxtVp/ + ]^2Dd''dl J /2(r, p, r, p', t)dp' = 0. 



(70) 
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which contains quasi-spin distribution function (129|) . where we have operators instead 
of the Pauh matrices a^, and, in the same way, we obtain kinetic equation for the spin- 
distribution function 

+ ^TsD VrVp j M^ir, p, r, p', t)dp' = 0. (71) 



where 



C°^(r, V,t) = \j (^V^* {R, t) [5{v - r„,)5(p - D„,) + 5(p - D„)5(r - r„)^ 



X 



xayiilj{R,t) + h.c.\dR, (72) 

and 



T,n = ^j dr{rf^. (73) 

T2D is the interaction constant for short-range interaction in two-dimensional system. We 
have obtained additional one particle function C"^, which unlike spin-1/2 case does not 
represents via / and S". As for electrically polarized particles, in section (V), appearing of 
additional function make truncation more difficult. We have suggested some approximation 



for hydrodynamical description (8|, j42| 



59|, where we have considered collective wave 



dispersion. We can suggest analogous approximations, but developing method allows us to 
perform truncation using approximate form of the wave function. 



IX. CONCLUSION 



Methods of studying of the quantum plasma have been developing last decade. The task 
of the interparticle interaction account reveals as a nontrivial problem. However we can 
conclude that this problem has been successfully solve in nonrelativistic in semi-relativistic 
quantum hydrodynamics. In classical physics kinetics method is more popular and profound 
method of the plasma studying. Thus it is worthwhile to derive the quantum kinetic theory, 
and it is also important to trace the interparticle interaction from microscopic level, as it 
was made in quantum hydrodynamics. First of all the kinetics allows to consider thermal 
effects. In most cases, when we have deal with the quantum plasma, quantum properties 
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might reveal at low temperatures. Thus using of the hydrodynamics in quantum case is 
more useful than in classical one. Nevertheless a some thermal spreading might leads to 
new effects in quantum plasma, and we have to have quantum kinetic theory to get correct 
description. 

This program has been realized in the paper. We have obtained quantum kinetic descrip- 
tion of different physical systems and presented corresponding kinetic equations in various 
approximations. 

We have suggested new definition for the quantum distribution function which gives us 
possibility to derive a chain of the kinetics equations directly from many-particle Schrodinger 
equation including explicitly inter-particle interaction. This derivation is the direct gener- 
alization of the many-particle quantum hydrodynamic method. This definition changes at 
account of the vector potential of the external electromagnetic field by evident explicit re- 
placement of the momentum operator p„ = —tKVn on D„ = — 2^V„ — e„ A„/c. In the 
quasi-classical limit it gives the Lorentz force in the usual form, and it also corresponds to 
the quantum hydrodynamics of spinning particles. Developed description corresponds to 
the Vlasov equation in the classical limit. 

We have derived kinetic equations for charged spinning particles (quantum plasma), 
neutral particles with short-range interaction, electrons in graphene, electrically polarized 
molecules, and graphene excitons. We have shown that in the system of spinning particles we 
have to use set of two kinetic equations, one for usual distribution function and the another 
for spin-distribution function which arise in the equation of distribution function evolution. 
We have obtained closed set of kinetic equation which allow to use them for studying of 
particular problems. 

Appendix A 

Here we will proof following formula ([8]), which we rewrite here 



we will make including presence of an external magnetic field, that corresponds to the fact 
that we should use general definition of the distribution function (123|) . Distribution function 
/ defined by formula fl23l) . or in the absence of the magnetic field and the rotational electric 




(Al) 



24 



field the distribution function lias more simple form ([5]). Hydro dynamical current j(r,t) has 
form Q, 3, B, B, Isl 



(A2) 



Now we are ready to demonstration of the formula ([8]) or flASp . 

p/(r,p,t)c/p = 

= \jpj $^||v'*(^'5(r-r„)5(p-p„) + 5(p-p„)5(r-r„)^^/^ + c.c.jdp 

= i / p / E (^(^ - ^")^(p - P-) + '^(p - p-)'^(^ - ^«)) ^ 

+ ^5(r - r„)5(p - p,J + 5(p - p„)5(r - r„)^ ■ dp 

j^jy^* (^^(r-r„)(5(p-p„) + (5(p-pJ5(r-r„)^pV' 

+ ^5(r-r„)5(p-p„) + (5(p-p„)5(r-r„)^(p^)* ■?/'jrfp 

Integrating on momentum p we replace p on p^ due to the delta functions 5(p — p„) and 
find ^ 

j p/(r,p,t)dp = ^ j |^lUp„^)*(^5(r-r„) + 5(r-r„))^ 



y{Y - r„) + 5(r - r„) j p„V^ j dp = j(r, t) 
Analogously, in the presence of the external magnetic field, we can show that 

= / p/(r,p,t)cip. 



(A3) 



where, from hydro dynamic point of view 

j^^' ^) = / E ^^(^ - ((DV^)*^^ + r^^^ dR. (A4) 
In formula ([9]) the hydrodynamic kinetic energy was mentioned, it's explicit form is 

C ^ 1 / 2 2 \ 

£(r, t)= Yl ^^(^ - ( (D V)> + V^*D V J rfi?. 

^ n=l ^ ^ 



(A5) 
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Formula ([9]) can be proved analogously to the proving of formula ([8]) or (1A3|) . 

The hydrodynamical spin-current appears in the formula fH2|) . where the second term we 
can rewrite as 

^ j d^ip^S^ir, p, t))dp = j /^"(r, p, t)dp = d^r\v, t), (A6) 



where 



(A7) 



t can be proved that formula (1A7I) corresponds to the QHD definition of the spin-current 



J, 14, [601. 



1. discussion of some methods derivation of kinetic equations 

In nineteen century, Boltzmann suggested his famous kinetic equation to describe an 
evolution of a gas of interacting particles. To describe interaction he introduced the collision 
integral. Since then the idea of using collision integral become widely used. It has been used 
even for plasma description, along with self-consistent field. Fundamental meaning of the 
collision integral we can understand by reading famous Landau's course of theoretical physics 
6ll |. In last decade a Boltzmann type equation has been used to description of ultracold 
quantum gases to describe a part "existing in excited states" that corresponds to the part 



of particles at non-zero temperature 62]- 65|]. A derivation of kinetic equation from NLSE 



was made in Ref. 



66| using Wigner distribution function. 



Let's pay attention to the ideas presented in Landau's course of theoretical physics |61 |. 
At description of non-equilibrium states of a Fermi- or Bose-liquid one uses quasi-particle 
distribution function n. This function depends on momentum, coordinate and time, and 
satisfies a transport (kinetic) equation. One takes the transport equation in the form of 

= /(n), (A8) 

where I{n) is the collision integral and d/dt is an operator which has following form: 

d d ^ d d 
dt dt dr ^(9p' 

where r and p should be find from Hamilton's equations. 
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Using of such representation to find an evident form of a transport equation is not cor- 
rect, because transport equation is specified in six dimensional space, where r, p and t are 
independent variables and quantities r and p are not specified and have no sense. 

Consideration of transport equation as a continuity equation in six dimensional phase 
space was very useful, it was a way of understanding of meaning of transport equation by 
many physicists. But kinetics chains have structure different from structure of hydrodynam- 
ics one. 

We have made criticism of the way of obtaining the left-hand side of the transport equa- 
tion, but one of the crucial points of this paper that we must not use form ( 1A8|) for transport 
(kinetic) equation. Moreover, we should not expect existence of collision integral in kinetic 
equation for description of short range interaction of neutral particles. An example of kinetic 
equation for neutral particles without collision term is presented in Ref. 66|. 



New kinetic equation gives us another method for study of particle dynamics. This equa- 
tion is derived for a system of neutral interacting particles. One does not contain collisions 
integral. It differs from Boltzmann equation. Equation contains interaction via term con- 
taining the trace of two-particle concentration /2 and interaction constant. The equation 
contains field variables only. In this case we can call it field kinetic equation. We obtain 
equation for the case of particles with short-range interaction. We do not use conditions 
of that interaction is weak or concentration is small. We found equation as for Bose as for 
Fermi particles, difference between Bose and Fermi particles reveals at further calculations 
of /2. It is interesting that we do not need to know cross-section for particles collisions as it 
is in Boltzmann equation. In this way we are free to calculate various properties of particles 
system without direct knowledge inter-particle potential of interaction. Equation ( 157|) con- 
tains information about inter-particle potential of interaction via interaction constant. The 
last one is the integral of inter-particle potential of interaction. The interaction constant 
can be found from measuring of dispersion dependence calculated from our equation. We 
do not need to know information about scattering length or cross-section before solving of 
the equation. 
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